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| Abstract 

^"J , In this paper, wc propose a stochastic model for a file-sharing pecr-to-peer network 

' which resembles the popular Bit Torrent system: large files are split into chunks and a 

peer can download or swap from another peer only one chunk at a time. We prove that 
the fluid limits of a scaled Markov model of this system are of the coagulation form, spe- 
cial cases of which are well-known epidemiological (SIR) models. In addition, Lyapunov 
stability and settling-time results are explored. We derive conditions under which the 
Bit Torrent incentives under consideration result in shorter mean file-acquisition times 
for peers compared to client-server (single chunk) systems. Finally, a diffusion approxi- 



, mation is given and some open questions are discussed. 
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> ■ 1 Introduction 

• rH 

X 

Peer-to-peer (p2p) activity continues to represent a very significant fraction of overall In- 
ternet traffic, 44% by one recent account [t]. Bit Torrent [H El El [HI [18l El HH] is a widely 
deployed p2p file-sharing network which has recently played a significant role in the net- 
work neutrality debate. Under Bit Torrent, peers join "swarms" (or "torrents") where each 
swarm corresponds to a specific data object (file). The process of finding the peers in a 
given swarm to connect to is typically facilitated through a centralised "tracker" . Recently, 
a trackerless Bit Torrent client has been introduced that uses distributed hashing for query 
resolution [16] , 

For file sharing, a peer is typically uploads upload pieces ( "chunks" ) of the file to other 
peers in the swarm while downloading his /her missing chunks from them. This chunk swap- 
ping constitutes a transaction-by-transaction incentive for peers to cooperate (i.e., trading 
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rather than simply download) to disseminate data objects. Large files may be segmented 
into several hundred chunks, all of which the peers of the corresponding warm must collect, 
and in the process disseminate their own chunks before they can reconstitute the desired 
file and possibly leave the file's swarm. 

In addition to the framework in which data objects are segmented into chunks to promote 
cooperation through swapping, there is a system whereby the rate at which chunks are 
uploaded is assessed for any given transaction, and peers that allocate inadequate bandwidth 
for uploading may be "choked" [14} 117] . Choking may also be applied to peers who, by 
employing multiple identities (sybils), abuse Bit Torrent's system of allowing newly arrived 
peers to a swarm to just download a few chunks (as they clearly cannot trade what they 
simply do not as yet possess). Bit Torrent can also rehabilitate peers by (optimistically) 
unchoking them. In the following, we do not directly consider upload bandwidth and related 
choking issues. 

In this paper, we motivate a deterministic epidemiological model of file dissemination for 
peer-to-peer file-sharing networks that employ Bit Torrent-like incentives, a generalisation of 
that given in jlojE Our model is different from those explored in |15 ]I211[T8] for BitTorrent, 
and we compute different quantities of interest. Our epidemiological framework, similar to 
that we used for the spread of multi-stage worms [12] . could also be adapted for network 
coding systems. In [9], the authors propose a "fluid" model of a single torrent /swarm (as 
we do in the following) and fit it to (transient) data drawn from aggregate swarms. The 
connection to branching process models [21] [8] is simply that ours only tracks the number 
of active peers who possess or demand the file under consideration, i.e., a single swarm. 
Though our model is significantly simpler than that of prior work, it is derived directly 
from an intuitive transaction-by-transaction Markov process modelling file-dissemination 
of the p2p network and its numerical solutions clearly demonstrate the effectiveness of 
the aforementioned incentives. A basic assumption in the following is that peers do not 
distribute bogus files (or file chunks) |20j. 

This paper is organised as follows. The Markovian model is developed in detail in Sec- 
tion [2j A proof of its fluid limit is described in Section [3] (and in the Appendix) including 
Lyapunov stability and settling time results. The behaviour of the limiting ODE is stud- 
ied in Section O for specific examples. In Section [5j we derive conditions under which the 
Bit Torrent-like system has improved performance (smaller mean time to completely acquire 
the file by a peer) compared to a system of pure client-server (no chunk-swapping) inter- 
actions. A diffusion approximation is given in [6J The paper concludes with a discussion of 
open problems in Section [7] 

2 The stochastic model 

We fix a set F (a file) which is partitioned into n (on the order of hundreds) pieces called 
chunks. Consider a large networked "swarm" of N nodes called peers. Each peer possesses 
a certain (possibly empty) subset A of F. As time goes by, this peer interacts with other 
peers, the goal being to enlarge his set A until, eventually, the peer manages to collect all 
n chunks of F. The interaction between peers can either be a download or a swap; in both 
cases, chunks are being copied from peer to peer and are assumed never lost. A peer will 
stay in the network as long as he does not possess all chunks. After collecting everything, 
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sooner or later a peer departs or switches off. By splitting the desired file into many chunks 
we give incentives to the peers to remain active in the swarm for long time during which 
other peers will take advantage of their possessions. 

2.1 Possible interactions 

We here describe how two peers, labelled A, B, interact. The following types of interactions 
are possible: 

1. Download: Peer A downloads a chunk i from B. This is possible only if A is a strict 
subset of B. If i E B then, after the downloading A becomes A' = A U {i} and but B 
remains B because it since it gains nothing from A. Denote this interaction as: 

(A <- B) ~» (A',B) 

The symbol on the left is supposed to show the type of interaction and the labels 
before it, while the symbol on the right shows the labels after the interaction. 

2. Swap: Peer A swaps with peer B. In other words, A gets a chunk j from B and 
B gets a chunk i from A. It is required that j is not an element of A and i not an 
element of B. We denote this interaction by 

(At5fl)v» (A',B') 

where A' = A U {j}, B' = B U {i}. We thus need A\B ^ and B \ A ^ 0. 

2.2 Notation 

The set of all combinations of n chunks, which partition F, is denoted by &(F), where 
\&{F)\ = 2™ and the empty set is included. We write A C B (respectively, AC- B) when A 
is a subset (respectively, strict subset) of B. We (unconventionally) write 

An A' when Ac A' and \A' - A\ = 1. 

If A H B = 0, we use A + B instead of ^4 U -B; if 5 = {6} is a singleton, we often write A + b 
instead of A + {&}. If A C B we use B — A instead of B \ A. We say that 

A relates to B (and write A ~ B) when A C B or B C A; 

if this is not the case, we write A ^ B. Note that A </< B if and only if two peers labelled 
A, B can swap chunks. The space of functions (vectors) from &{F) into Z + is denoted by 
"L^ F \ The stochastic model will take values in this space. The deterministic model will 
evolve in R + . We let E Z + be the vector with coordinates 

e%:=l(A = B), B E &>(F). 

For x E or (F) we let |x| := Ea6^(F) 1^1- If ^ C ^(F) then the j^-face Kf of 

R+ {F) is defined by Hf := {x E rJ (F) : Ea £ ^ ^ A = °> IIb^ ^ B > 0}. 
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2.3 Defining the rates of individual interactions 



We follow the logic of stochastic modelling of chemical reactions or epidemics and assume 
that the chance of a particular interaction occurring in a short interval of time is proportional 
to the number of ways of selecting the peers needed for this interaction [13]. Accordingly, 
the interaction rates must be given by the formulae described below. 

Consider first finding the rate of a download A <— B, where A C B, when the state of 
the system is x G zf* . There are x A peers labelled A and x B labelled B. We can choose 
them in x A x B ways. Thus the rate of a download A <— B that results into A getting some 
chunk from B should be proportional to x A x B . However, we are interested in the rate of 
the specific interaction (A <— B) (A', B), that turns A into a specific set A' differing from 
A by one single chunk (A C A'); there are \B — A\ chunks that A can download from B; the 
chance that picking one of them is 1/\B — A\. Thus we have: 

{x B 
the rate of the download (A <— B) (A' , B) equals X Tg — 
as long as A C A 1 C B, 

where j3 > 0. 

Consider next a swap A B and assume the state is x. Picking two peers labelled A 
and B (provided that A B) from the population is done in x A x B ways. Thus the rate of 
a swap A ^ B is proportional to x A x B . So if we fix two chunks i G A\ B,j G B \ A and 
specify that A' = A + j, £?' = B + z, then the chance of picking i from A \ B and j from 
£ \ A is 1/|,4 \B\\B\ A\. Thus, 

{X^X^ 
the rate of the swap (A^ B) (A' , B') equals ^ rjXjjjTjfT^p 
a long as An A', B C B', A' - A C B, B' - B C A, 

where 7 > 0. 

2.4 Deriving the Markov chain rates 



Having defined the rates of each individual interaction we can easily define rates q(x, y) of 
a Markov chain in continuous time and state space Z + v ' as follows. 

Define functions \a,A' ^A.B '■ 



x c 



E \C-A\ 
c-.cdA' 1 1 



1(A C A') (la) 



4 R 

X X 



Consider also constants 5 > and a" 4 > for A G ^(i 7 ), i.e., a G . The transition 
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rates of the closed conservative Markov chain are given by: 



Xa,A'(x), ify = x-e A + e A > 

y = x - e A - e B + e A > + e B > 
An A',B c B',A' - Ac B,B' - B c A, 



HA,b{x), if 



a it y = x + e A 

5x F \i y = x — ep 

0, for any other value of y 7^ x, 

where x ranges in zf*^. 

A little justification of the first two cases is needed: that g(x, x — e A — e B + e A i + 
e_B') = /m,_bOe) is straightforward. It corresponds to a swap, which is only possible when 
A \Z A',B \Z B',A' - A C B,B' - B C A. The swap rate was defined by (SR). To see 
that q(x,x — e A + e A >) = Xa,A'( x ) we observe that a peer labelled A can change its label 
to A' Z\ A by downloading a chunk from some set C that contains A' , so we sum the rates 
(DR) over all these possible individual interactions to obtain the first line in ([2]). We can 
think of having Poisson process of arrivals of new peers at rate \a\, and that each arriving 
peer is labelled A with probability o^/|a|. Peers can depart, by definition, only when they 
are labelled F and it takes an exponentially distributed amount of time (with mean 1/5) 
for a departure to occur. Thus, q(x,x — ejr) = Sx F . We shall let Q denote the generator of 
the chain, i.e. Qf(x) = ^2 y (f(y) — f(x))q(x,y), when / is an appropriate functional of the 
state space. 

Definition 1 (BITTORRENT [xq, n, a, (3, 7, 51 ). Given xq G zf**^ (initial configuration), n = 

\F\ G N (number of chunks), a G Jjf* (arrival rates), f3 > (download rate), 7 > 
(swap rate), 5 > (departure rate) we let BITTORRENT [xo, n, a, /?, 7, <5] 6e a Markov chain 
(Xt,t > 0) w'i/i transition rates ([2]) and Xo = xo- VFe say i/iai i/ie chain (network) is open 
if a A > for at least one A and 5 > 0; it is closed if a A = for all A; it is conservative if 
it is closed and 5 = 0; it is dissipative if it is closed and 5 > 0. 

In a conservative network, we have q(x,y) = if \y\ 7^ |x| and so \Xt\ = \Xq\ for all 
t > 0. Here, the actual state space is the simplex 

{x G zf (F) : |x| = N}, 

where N = \Xq\. It is easy to see that the state ep is reachable from any other state, but all 
rates out of ep are zero. Hence a conservative network has ep as a single absorbing state. 

In a dissipative network, we have \X^\ < \Xq\ for all t > 0. Here the state space is 

{x G Zf {F) : |x| < iV}, 

where iV = \Xq\. It can be seen that a dissipative network has many absorbing points. 

In an open network, there are no absorbing points. On the other hand, one may wonder 
if certain components can escape to infinity. This is not the case: 

Lemma 1. If a F > then the open BITTORRENT [x, n, j3, 7, a, 51 is positive recurrent Markov 
chain. 
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Proof, (sketch) If a F > 0, 5 > the Markov chain is irreducible. The remainder of the 
proof is based on a the construction of a simple Lyapunov function: 

V(x) := \x\, 

for which it can be shown that there is a bounded set of states K such that 

sup(Qy)(x) < 0. 

Perhaps the easiest way to see this is by appealing to the stability of the corresponding ODE 
limit; see Theorem Q] below and [7]. Q 

2.5 Example: n = 1 

Let us take the special case where the file consists of a single chunk (n = 1). The state here 
is x = (x ,x l := x F ). The rates are: 

q{(x ,x l ),{x + l,x 1 )) = a 

c/((x ,x 1 ),(xV 1 + l)) =a l 
q({x , x 1 ), (x - 1, x 1 + 1)) = Px x 1 

g((x ,x 1 ),(x ,x 1 - 1)) =5x\ (3) 

If a = q 1 = 0, this is the stochastic version of the classical (closed) Kermack-McKendrick 
(or susceptible-infective-removed (SIR)) model for a simple epidemic process [3]. Its absorb- 
ing points are states of the form (x , 0). In epidemiological terminology, x 1 is the number of 
infected individuals, whereas x is the number of susceptible ones. Contrary to the epidemi- 
ological interpretation, infection is desirable, for infection is tantamount to downloading the 
file. 



3 Macroscopic description: fluid limit 

Analysing the Markov chain in its original form is complicated. We thus resort to a first- 
order approximation by an ordinary differential equation (ODE). 

Let v(x) be the vector field on R^T^ with components v A (x) defined by 



Or) = a A -x A ((3<pj{x)+^f(x)) 

ip A (x)x B sr~^ ip^' B (x)x B 

_ + \B \ A\ ^ 2.*/ 

B-.ACB 1 V 1 B:A<£B 



^ ^ 1 + \B\A\ ' 1 + \B\A\ K h y ' 



where 



^(z):=J> B , ^(*):=J>* 



bda b^a 



a<=A a£A(lB 
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Consider the differential equation 

x = v(x) with initial condition xq. (6) 

Consider the sequence of stochastic models BITTORRENT [X N> o, n, Nat, jj, jj, 51 for N G N 
and let X N j be the corresponding jump Markov chain. 

Theorem 1. There is a has a unique smooth (analytic) solution to {6|), denoted by xt for 
t > 0. Also, if there is an xq £ K + such that X Nj q/N — ► xq as N —> oo, then for any 
T, e > 0, 

lim P( sup {N^Xxj - x t \ > e) = 0. 

N^oo 0<t<T 

Proof. Let M be the set of vectors — e^, e^, — + e^', — — ee + &A' + ^b', where 
A,B e ^(F) and A C A', 5 C 5'. From ©, we have that <?(x,y) = if y - x AT. 
Introduce, for each £ £ A/", a unit rate Poisson process 3>£ on the real line, and assume 
that these Poisson processes are independent. Consider the Markov chain {Xt) for the 
BITTORRENT LXq, n, a, (3, 7, 51 . Its rates are of the form 

q(x,x + () =Q ( (x), CGAA, (7) 

where Qc_(x) is a polynomial in 2 n variables of degree 2, and which can be read directly from 
([2]); its coefficients depend on the parameters a, (3, 7, 5. We can represent [131 [5] (Xt) as: 



Ii=X +V($c( [ Q c (X s )ds). 



Consider now the Markov chain -kX N j corresponding to to BITTORRENT [X N) q, n, (Na A ), /3/N,j/N, 51 . 
The transition rates for jjX N j are 

q(x/N,(x + ()/N) =NQ c (x/N), x E zf {F) , ( £ M, 

and 0, otherwise. Here, Q^(x) is the polynomial defined through ([7]) and ([2]) and we now 
assume that its variables range over the reals. Therefore, ij^j can be represented as 

±x N , t = 1a^ + £ C>c (n f Qd^x»,s)ds) . 

Define xt by the (deterministic) integral equation 

Xt = Xq + ^ C / Q((x s )ds (8) 

and assume that it is unique for all £ > 0. Fix a time horizon T > and let 

B := max \xA, 
t<T 1 h 

M c := max 

|a;|<B 

:= sup — — j p 

|x|,|j/|<B F ~~ 2/1 

t n := inf{t > : \X N>t \ > NB}. 
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We then have: 



A N , t := ^fr-xt = ^r- x ° + J2<\jj®c( N [ Qd X »,s/N)ds\ - f Q C (X N , S /N)ds 

£ejV L V Jo J Jo 

+ E C / (Qd^,./N) - Q c (x s ))ds 

TTfr Jo 



Suppose that t < T A t n . Then, for all s <t, 

\Q C (X N>S /N)-Q c {x s )\ <L C |A 

So, if we let 



Xjv,q 
?iv,t •— — ^ ^0 



+ *<( N J QdX^Md*) - N J* Q C {X N , S /N)ds 



we have, by the Gronwall-Bellman lemma, that 



Let 



IAjv.il < |^ it |exp(t^ |<|L C ), Xt<TAT N . 



&Ut) := sup|$ c (s) - s\. 

s<t 



We recall that, as N — > oo, 

^((iVi) - a.s. (9) 

If s < t < t n , we have X N . S /N < B (definition of r N ) and so Q^(X NtS /N) < M^, implying 
that 

sup \g Nit \<\^- X0 \ + Y,\C\^l(NM c T) 

t<TAr N ^ JV 

which converges to zero, a.s., due to ([9]). Since 

sup \A N)t \ < sup \£ N ,t\ exp (T ICI^c)) 

t<TAr]v t<TArjv fgAr 



we have 



Now observe that 



sup |A N ,t| — > 0, a.s. 

t<T/\T N 



P(t n <T)<P( sup \X N)t \ > NB) 

t<TAr N 

<P( sup |Ajv. t | + sup \x t \>B)—>0. 

t<TAr N t<TAr N 

So we have sup t<T |Ajv,i| — ► a.s. 

To show that Xt, defined via (JSj) , satisfies the ODE x = v(x) with v given by is a 
matter of straightforward (but tedious) algebra, see the Appendix. 
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Uniqueness and analyticity of the solution of the ODE is immediate from the form of 
the vector field (its components are polynomials of degree 2 and hence locally Lipschitzian). 

To show that the trajectories do not explode, we consider the function 

V(x) :=£V. 
A 

It is a matter of algebra to check that 

(VV(x),v(x)} = £V(x) = J2 aA ~ 5xF 
A A 

which (since 5 > 0) is negative and bounded away from zero for x outside a bounded set of 
R + containing the origin. We then apply the Lyapunov criterion for ODEs to conclude 
that xt is defined for all t > and this justifies the fact that we could choose an arbitrary 
time horizon T earlier in the proof. □ 

Comment: The quantities defined in ([5]). have physical meanings as follows: 

no. of peers from which an ^4-peer can download, 

no. of peers an A-peer can swap with, 

= no. of peers which can assume label A after a download, 

a = no. of peers which can assume label A after a 5-peer swap. 

It is helpful to keep these in mind because they aid in writing down the various parts of 
v(x), again, see the Appendix. 



<rf(x) := £ x B 



bda 

B 



■= Yj x 

B^A 



x 

A- 



aeAnB 



4 Behaviour of the limiting ODE 

Concerning the ODE x = v(x) we consider again three cases, just as in the stochastic 
model: an open system (a A > for at least one A and S > 0), a closed dissipative system 
(a = for all A and S > 0), and a closed conservative system (a A = for all A and 

5 = 0). Qualitatively, the behaviour is different in each case. In this section, we will try to 
exemplify this behaviour by means of examples. 

4.1 The ODE in the absence of BitTorrent incentives 

Absence of BitTorrent incentives means that the file is not split into chunks, i.e. n = 1. 
Thus, the Markov chain is Xt = (X t ,X^ := Xf), i.e. two-valued A £ {0, 1}. The rates for 
this case were reported earlier in ([3]). To find the fluid limit, we use ([5]) and (j4]), keeping in 
mind the interpretation of each of the terms in ((5]). 

For A = 0, we have: the number of peers from which an 0-peer can download from 
is ip d (x) = x 1 ; the number of peers that can swap with an 0-peer is ipf(x) = 0; since no 
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peer can assume value after an interaction, we have ip d (x) = ijj s ' (x) = 0. Hence, in the 
formula for (JH) for v (x) only the first parenthesis survives: 

v {x)=a - (3x x l . 

For A = 1 we have: f d {x) = fl(x) = 0, ^ l d {x) = x , i/jI ,0 (x) = 0. Thus, 

v 1 (x) = a 1 +[3x x 1 -5x\ 

Here, by simply letting x := x and y := x 1 , the ODE is 

x = a — j3xy 

y = a 1 + j3xy — by. 

If a = a 1 = 5 = (closed conservative system), we have x + y=constant, say = 1, and 

x = —(3x(l — x), 
the solution of which is the logistic function, 




12 3 4 



If a = a 1 = 0, but 5 > (closed dissipative system) we have the classical deterministic 
SIR epidemic 

x = —j3xy 

y = fay - $y 

the integral curves of which can be found by solving 

^ = — - l 

dx (3x 

whence it follows 



y = (x + y ) + - log(x/x ) - x: 
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Assume xq + yo < 1. Notice that the integral curves are monotonic as long as they start 
from a point (xo, yo) with xo > 5/(3. They eventually converge to a point of the form (x*,0) 
where 

(zo + Vo) + p log(x*/x ) - x* = 0. 

In other words, (xt,yt) — ► (x*,0) as t — > oo and, in fact, does so exponentially fast. 

Notice that if the initial state is in the interior of the positive orthant then the boundary 
cannot be reached in finite time. This is in contrast to the corresponding stochastic model 
which can reach the boundary in finite time with positive probability. In fact, in the open 
case, it will reach the boundary in finite time with probability 1 due to (positive) recurrence. 
This remark is generic: it applies to any dimension. 

Next, consider the open system case assuming, for simplicity, that a = A > 0, but 



a 



0. We have 



x = A — (3xy 
V = (3xy - 5y. 

Since 5 > 0, the Lyapunov function argument shows that the system is asymptotically 
stable. In fact, there is a unique asymptotic equilibrium which attracts all initial conditions 
(2:0,2/0) w ith yo > 0. This unique equilibrium is given by 

* 5 * A 
X = -^ V = T 

as can be seen by setting the right hand side of the ODE equal to zero [3]. It should be 
noticed that the trajectories can be spirals around (x*,y*). A typical situation is shown 
below. 



y i.o- 



\ \ \ \ \ 

\ \ \ \ \ 

\ \ \ \ \ 

\ \ \ \ \ 

\ \ \ \ \ 

\ \ \ \ \ 

\ \ \ \ \ 

\ \ \ \ \ 

\ \ N \ N 

X \ X x x 



In this vector field plot, we took (3 = 3, A = 5, 5 = 4, so (x*,y*) = (1.33,1.25). If yo = 
then the trajectory converges to infinity. Indeed, if no peers initially possess the file and no 
such peers ever show up, then the only thing that can happen is an accumulation, at rate 
A, of ever demanding peers. The remedy is, clearly, the imposition of a 1 > no matter how 
small; then (x*,y*) is globally attracting for all initial states in the closed positive orthant. 

In this last example, we find that the eigenvalues of the differential of the vector field at 
(x*,y*) are complex conjugates if and only if A/3 < 4<5 2 , and so this is the condition for the 
spiralling of trajectories. 
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4.2 The ODE for n = 2 chunks 



Here A can take 4 values: 0, {1}, {2}, {1,2} := F. We work out the expression for each 
component separately. Since an 0-peer can only download from <£ (x) = x 1 -\-x 2 +x 12 peerg^, 
but no peer can download from it or swap with it, we have <pf(x) = ip {x) = ipf{x) = 0, 
and so 

v ( x ) = a - 0x {x l +x 2 + x 12 ). 

A 1-peer can download from tp\(x) = x 12 peers; it can swap with <fl(x) = x 2 peers; the 
number of peers that take value 1 after a download is i^\{x) + x ; and, since no peer can 
take the label 1 after a swap, we have i/jI ,B (x) = V.B. We thus have 

v 1 (x)=a 1 -x\(3x 12 + 1 x 2 )+(3ip 1 d (x) £ 

B-.ieB ' \ \ s\ 

= a l - x\f3x 12 + 7 x 2 ) + (3x (x l + ix 12 ). 

The expression for v 2 (x) is symmetric to that of v 1 (x). A 12-peer can neither download 
from or swap with anyone, so ip^ 2 (x) = Lp l2 (x) = 0; there are ip^ 2 (x) = x 1 + x 2 peers which 
can take label 12 after a download; since only a 2-peer can assume label 12 after a swap, we 
have ijj\ 2,l {x) = x 2 ; and, for the same reason, ipl 2,2 (x) = x 1 . Thus, 

12/ n 12 , fl ,12m 12 , -fipl 2,l (x)x 1 7^ 2,2 (x)x 2 12 



We thus have, 



+ |{1}\{1,2}| 1 + |{2}\{1,2}| 
a 12 + Pix 1 + x 2 )x 12 + -y(x 2 x l + x 1 ^ 2 ) - 5x 12 . 



x = a -(3x (x 1 + x 2 + x 12 ) 

i 1 = a 1 - x l (Px 12 + 7 x 2 ) + (5x {x l + \x 12 ) 

± 2 = a 2 - x 2 {(5x 12 + 7X 1 ) + (3x {x 2 + \x 12 ) 

^_19 19 . m 1 . 9i 19 . „ 19 



^12 _ a 12 + ^^1 + x 2 )x X2 + 2lx X x 2 _ §x 12 



Case 1: closed conservative system. Consider the n = 2 case again, and assume that 
a = for all A, 5 = (a closed conservative system). Assume further that 7 = 0. Let 

x = x , u = x 1 + x 2 , to = x 12 . 

We see a reduction in dimension from 4 to 3 (owing to that fact that the sum of the 
coordinates is constant) and a further reduction from 3 to 2 (owing to the fact that, for 
7 = 0, the vector field depends on x , x 2 through their sum). We have: 

x = —f3x(u + w) 

u = —(3uw + (3x{u + w) 

w = (3uw. 

On assuming that 

X + u + w = 1 

^Where x 12 = x F and indexes 1 and 2 represent the chunks that form a partition of F. 
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and eliminating the variable w we obtain 



x = — /3x(l — x) 

ii = f3u 2 - 0u(l - x) + @x(l - x). 

We can qualitatively see the behaviour of this ODE by looking at the vector field in the 
x — u plane. A typical picture is as follows: 



~- v, v, V N N 
~ •>« WN 
^ ~- ^ w N 

— ~~ ^ ^ V N N 
~- ^ ^ x V N, V 
* ^ ^ V *x V 

— ^ ^ ^ X V. V 
^ ~v ^ W \ 

^ ^ -v. V. -x V 
^ ^ ^ X V V 
~- x ^ X v, v. V 
~- x X *x X V V 
^ X V X N V V 
■x X X V N V V 

V S N \ W\ 
N V X V V N N 

V N V S \ S \ 

wwww 

\ \ N \\ \\\ 



(').:■) 



0.4 



\\ 
\\ 

\\ 
V N. 

\S 
\\ 
\\ 



\\ 
\\ 
\\ 

0.5 



The first is an autonomous equation, encountered earlier. Its solution is x t = xo/(xq + (1 



•i'O 



3t 



. Letting 



/:=1 



w 



X + u 



we have 



f = Pf(l-f)+P(l-f) Xt 
Assume that wo > 0, i.e. fo < 1. The solution of the last ODE is: 

/o(l-/o) _1 + a?o^ 



Hence 



ft 



w t 



x + (1 - x Q )ef* + / (1 - /o)" 1 + x pf 
xq + (1 - x )e (St 



xo + (1 - x )ef 3t + x f3t + (1 - w )w 

We have lim^oo wt = 1, as expected. From this, we can estimate the time r = t(xq,wq,£) 
required for the system to reach the set {(x, u,w) : w > 1 — e}, where < e < 1 is a given 
(small) number. This time can be translated into an estimate for the expected time required 
for the original stochastic system to be absorbed by the state where all peers possess the 
full file. We shall not attempt the justification of this statement here. The time r is the 
solution of the transcendental equation 



1 - Wq 



+ XqPt 



: (x + (l-x )e /3T ). 



w 1 
Typical behaviour of this time as a function of xq is as follows: 



13 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



The three graphs correspond to varying values of (3 ranging from 1 (top curve) to 5 (bottom 
curve) and to e = 0.001, wq = 0.1. 

Case 2: closed dissipative system. Consider again the n = 2 case with no arrivals, 
with swap rate 7 equal to zero, but with departure rate 5 > 0. Using the same notation as 
above, we have 

x = —j3x{u + vo ) 

u = —(3uw + (3x(u + w) 

w = (3uw — 6w. 

To substitute out one parameter, change the time variable to s = fit, let p = 5/(3, and write 
x' for dx/ds, so that 

x 1 = —x(u + w) 

v! = —uw + x(u + w) 

w = uw — pw. 

Assume 

x + u + w = 1, 

so xt + ut + wt < 1 for all t > 0. We cannot eliminate the variable w now since there is no 
obvious conserved quantity, but we can study the equilibria of the system. It seems that the 
ODE has no obvious analytical solution. However it can either be integrated numerically or 
its solution in terms of power series can be easily found. Setting the right hand side equal to 
zero we see that, necessarily, w = 0, which leaves us with the possibility xu = 0. So points 
of the form 

(x, u, w) = (x, 0, 0), (0, u, 0) 

are equilibria. But not all of them are stable. For example, any point of the form (0, u, 0) 
with u > p is an unstable equilibrium. Specifically, the differential of the vector field is given 
by 

(—(u + w) —x —x \ 
u + w x — w x — u\ 
w u — p) 
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/O \ 

Evaluated at (0, u, 0), it gives Dv = — u , and this has eigenvalues 0, —u, u — p. 

\0 u-p/ 

(0 —x — x\ 
x x I , and this has eigenvalues 0, x, — p. 
p) 

Therefore, the only stable equilibria are points of the form 

(x, u, w) = (0, u, 0), < u < p. 
In terms of the original variables, the stable equilibria are 

(x ,x l ,x 2 ,x 12 ) = (0,x\x 2 ,0), < x 1 +x 2 < p. 

This is as expected: since there is no swapping (7 = 0), the system eventually settles to a 
situation where there are peers with label 1 and peers with label 2. Had 7 been positive, 
x 1 , x 2 could not have simultaneously been positive in equilibrium. 



Case 3: Open system. Consider the situation as in Case 2, but add arrivals of peers 
(known as seeds) possessing the full file. Choosing variables appropriately, we have 

x' = —x(u + w) 

v! = —uw + x(u + w) 

w' = A + uw — pw. 

In terms of the original variables, we here have a = a 1 = a 2 = 0, a 12 = X(3 > 0. Here we 
see that the system eventually settles to the stable equilibrium 

(x,u,w) = (0,0, X/p). 

We can easily see that the eigenvalues of the differential of the vector field at the stable 
equilibrium are both real and negative: —X/p and — p (the first one has algebraic multiplicity 
1 but geometric multiplicity 2). Hence there is no possibility of spiralling. 

w 
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4.3 Time to settle 

Estimating the time for the system to reach an equilibrium requires cooking up an appro- 
priate Lyapunov function. The obvious Lyapunov function used earlier gives a crude lower 
bound. 

Consider a general deterministic system (open or closed), i.e. the differential equation 
x = v(x) with v(x) given by (jlj). Let V(x) = \x\ = ^2acf xA - Then 

V = \a\ — 5x F ', 

where \a\ = ^2^cF a A is the total arrival rate. Using the naive inequality x F < V we obtain 
the following: 

Corollary 1. If the system starts from xq and if xt — ► x* then the time r r required for the 
trajectory to reach an r-ball centred at x* satisfies 

1 \x \ + \a\ 

T r > - log 



S \x*\ + r 

In the other direction, consider the last component v F (x) of the vector field. Since no 
seed (F-peer) can download or swap, we have <p F (x) = <pf(x) = 0. So 



v F 



(x) = a F + 01$ {x)x F + 7 ^ B {x)x B - 5x F 

BcF 

= ^ + 0^ x F x F ~i + 7 E E ^ /_! " 6xF 

jeF B<ZF ieB 

=: a F + v F (x) — 5x F , 

where v F (x) is a quadratic form with positive coefficients. 

Assume that the system is closed, so that a F = in particular. Let v F be an upper 
bound on v F (x), i.e., a (good) bound that depends on the initial state xq and the parameters 
P, 7- 

Corollary 2. If the system is closed and v F < 5, then 

1 x F 

T r < = log-2-. 

v F -5 r 

We conjecture that v F (x t ) decreases along the trajectory xt as long as v f (xq) < 8. If 

p 

so, the bound above is 1 , — F log — . 

v^_(xo)—o ° r 

5 An example of the evaluation of performance improvement 
in presence of Bit Torrent incentives 

We address the following question: When is it advantageous to split a file into chunks? 
In other words, assuming we fix certain system parameters (e.g., arrival rates), will peers 
acquire the file faster if the file is split into chunks? We attempt here to answer the question 
in a simple case only by using the deterministic approximation. Let A be the total peer 
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arrival rate. Let (5 be the download rate. Assume that only peers arrive exogenously. In 
the absence of Bit Torrent incentives, we have the single-chunk case 

x = X - (ix x l 
x 1 = (3x x 1 -Sx 1 . 

The globally attracting stable equilibrium is given by 

x* = {5/(3, X/6). 

Consider splitting into n = 2 chunks. Let x be the state of the system. Suppose that the 
new parameters are A = A, 6 = 5, (3, 7. Then 

= X -Px (x 1 +x 2 + x 12 ) 
x 1 = -x l (Px 12 + jx 2 ) + fix {x l + \x 12 ) 
x = -x 2 (Px 12 + 7s 1 ) + fix (x 2 + \x 12 ) 
i 12 = p(x l + x 2 )x 12 + 27^^ - 5x 12 . 

The new equilibrium is easily found to be 

u u A 
,/3 v A~ ' ~> ' 2' 2' 7 

where it is the positive number which solves 

q{u) = 0, 

and 

2 2/3A 2A , in . 
70 7 

To see this, set the vector field equal to zero and solve for x A , A 6 {0, 1, 2, 12 := i* 1 }. The 
simplest way to obtain the solution is by first adding the equations; this gives 

A - 6x 12 = 0, 

whence x 12 = \/5. Then add the middle two equations after setting u = x 1 + x 2 : 

-fiux 12 - 2jx l x 2 + f3x (u + x 12 ) = 0. 

Replace x 12 by X/6 and observe that, due to symmetry, x 1 = x 2 = u/2. This gives the 
quadratic equation q{u) = with q defined by (jlOp . Finally, the first equation becomes 

X-Px (u + X/5) =0, 

which is solved for x giving: 

X, A, 1 5 ,5 5 
x* = ^(u+ -) 1 = ^(-u+ 1) < ^. 

Thus: 
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Corollary 3. If > then 



x* < x 



*0 



So, by introducing splitting into chunks, we have fewer peers who have no parts of the 
file at all. Using Little's theorem (see below), this can be translated into smaller waiting 
time from the time a peer arrives until he gets his first chunk. 

Suppose now we are interested in determining how long it will take for a newly arrived 



peer to acquire the full file. On the average, a peer spends time equal to A 



before it 



exits the system. During last part of his sojourn interval (which is a random variable with 
mean 1/5), the peer possess the full file. It thus takes on the average A -1 |a;*| — 5" 1 for a 
peer to acquire the full file. Since we assume that A = A, 5 = 5, it suffices to show that 



> x* 



But 



x 



6_ A 
~0 + 5 



5 ,6 



-u 



0^ 



o 6 ,6 s-i 



A 

U ~~5 



recalling u > solves q(u) = 0. So, \x* \ - 

> q(u) := u 2 



f 5__5_ l_ 

^ f ^ K 0\ 



A 



1 )u- -u' 



\x*\ > if and only if 



,5 A, ,A As 



(") 



Define u as the unique positive number which satisfies 

q(u) = 0. 

Corollary 4. If > j3, a necessary and sufficient condition for \x*\ > \x*\ is u <u. 

This gives a set of non- vacuous conditions for achieving improvement of performance by 
the introduction of BitTorrent incentives!! It can be proved that if the parameters 0, S, 0, 
7 are fixed and if > then there exists a value Ao such that for all A < Ao the inequality 
u <u holds. To prove this, we observe that 



u < u 



q{u) > 



and study the behaviour of q(u) as a function of A in a neighbourhood of zero. 

We conjecture that an algebraic condition involving quadratics like q and q is valid for 
larger values of n also. 

To justify the use of deterministic approximation for estimating performance measures, 
and, specifically, the use of mean values, we need to show that as N — > oo, we can approx- 
imate stationary averages in the original stochastic network by equilibria of the resulting 



*The inequality conditions were mistakenly reversed in the corresponding results of [11] . 
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ODE. It is easy to show that the a.s. convergence to the ODE limit can be translated into 
convergence of the means, using a uniform integrability argument. Namely, 

—EX Nt ► x t > x*, 

iV N — >oo t— >oo 

where the second limit concerns the behaviour of the ODE alone. On the other hand, if we 
fix N and look at the asymptotic behaviour of the process j^X N: t as t — > oo, we have 

1 1 ~ 

where the law of Xjq is the stationary distribution of the chain (jjX N j)t>o- It can be proved 
that jjEX n N — > oo. Arguments for this will be considered in future work. More 

detailed estimates on the discrepancy between the stochastic and deterministic systems can 
be found in the recent survey paper [6|. 

We can also explain the use of \x*\ as a measure of the sojourn time in the system of 
a peer, by first using the approximation outlined above and then appealing to Little's law. 
This is as follows. 

Consider an open BITTORRENT IXq, n, a, /5, 7, 51 , i.e. \a\ > 0, 5 > 0. We know that the 
Markov chain (Xt) is positive recurrent and has thus a unique stationary distribution. It 
makes sense to assess the performance of the network by looking at steady-state performance 
measures, such as the mean time it takes for an 0-peer to become an F-peer (a seed). 
Consider then the process (Xt,t G M) defined to be a stationary Markov process with time 
index M and transition rates as those of (Xt). The law of the process (Xt,t G K) is unique. 
Let T A , k G Z be the times at which A-peers arrive (and, say, T A < < T^ 4 , by convention). 
These are the points of a stationary Poisson process in K with rate a A . Let W A be the 
sojourn time in the system of a peer arriving at time T A . Since, by assumption, a peer 
departs only after it has acquired the full set, the time W A is the sum of the times it takes 
for the peer to become a seed plus the time that the peer hangs out in the system after 
becoming a seed (the latter is an exponential time with mean 1/5). Clearly then, for all 
t G R, 

^X^ = ^HT k A <t<T k A + W k A ). 

bdA fcez 

Using Campbell's formula, we obtain 

EXo = a A E A W A , (12) 

where E A is expectation with respect to P^-the Palm probability of P with respect to the 
point process (T A , k G Z). 

In particular, with A = 0, and A = a , we have that 

E W O = jE\X \, 

which can be read as: the mean sojourn time of a 0-peer is, in steady state, equal to the 
mean number of peers in the system divided by the rate of arrivals of 0-peers. If TV" is a 
parameter of the process as in Theorem [1] then, A being proportional to N, we have that 
the right hand side converges to something that is proportional to | X I j clS required. 
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6 Diffusion approximation 



Using the functional central limit theorem for Poisson processes, we can prove, by standard 
methods, the following: Again consider the sequence BITTORRENTLX^Ojre, Na, 4?, jt,5] for 
N £ N, and let X N t be the corresponding jump Markov chain. Let (xt,t > 0) be the 
solution to the ODE x = v(x) with initial condition xq. Let 

Y N>t := VN(X Ntt /N - xt). 

Let W^, C £ -A/", be i.i.d. standard Brownian motions in R. Finally, define the (time- 
inhomogeneous) Gaussian diffusion process Y by 

dY t = C^Qd x t)dW ct + Dv(x t )Y t dt, 
where Dv(x) is the matrix of partial derivatives of v(x). 

Theorem 2. If \/N(X Ni t/N — xq) — >• as N — > oo, w/iere xo G R+ , i/ten i/ie kw ofY N 
(as a sequence of probability measures in D[0, oo) mi/i t/ie topology of uniform convergence 
on compacta) converges weakly to the law ofY. 

The proof of this theorem is omitted but the reader is referred to [13] for the relevant 
arguments. 

7 Final remarks, open problems and future work 

7.1 Rates of convergence 

We can obtain a computable rate of convergence of the stochastic model to the ODE by 
using a combination of large deviations techniques with the solution of two optimisation 
problems. The idea is basically implicit in the proof of Theorem [1] and this is the reason we 
wrote the proof explicitly in terms of the driving Poisson processes 

The first problem is so that we obtain an estimate of the maximum value of of Q((x). 
In the case of a closed network, this is a quadratic optimisation problem over the polyhedron 

{M < !}• 

The second optimisation problem is for an estimate for L^, which can be translated to 
an estimate for the norm of the gradient VQ^(x). In the case of a closed network, we can 
estimate this by solving a (large) number of linear programming problems. Of course, only 
estimates are needed. 

7.2 Conjectures 

The first one concerns the behaviour of v+(xt), and was stated at the end of Section [H The 
second concerning generalisation of Corollary U was stated in Section [5j 

The third conjecture is more vague: it basically says that we can evaluate the perfor- 
mance improvement by using a large number of chunks (say 100), by solving a number of 
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quadratic inequalities. To this end, it should be remarked that in a deterministic open net- 
work, the unique equilibrium x* can be found by solving n + 1 equations in re + 1 unknowns, 
provided that the rate of arrivals of ^4-peers depends on A through its cardinality alone: 

a A = Xk if \A\ = h chunks.. 

Indeed, by symmetry of the vector field, we see that 

x* A = x* B if \A\ = \B\. 

So if we define 

:= J2 x* A 

\A\=k 

we will have 

-l 

z k , if \A\ = k. 
Hence if we let 

V k {z°,...,z n ) := v A {x*), 

\A\=k 

where z and x* are related as above, then the equations we need to solve are 

V k (z°,...,z n ) = 0, k = 0,...,n. 

7.3 A reduction of dimension for a balanced ODE 

If we are interested not only in the equilibria but also in a more detailed study of the 
transient behaviour of the ODE, then we can obtain a rough idea (and bounds) by making 
the assumption of full symmetry, i.e., we assume that the arrival rates a A and the initial 
states x A depend on A only through \A\. Then the trajectory itself xf depends only on 
the cardinality of A and so we can reduce the ODE to an (re + l)-dimensional one. Such 
a symmetrised ODE can yield more detailed information on the time to reach a small 
neighbourhood of the equilibrium point. (Note that our bounds in §4.31 are very crude.) 




7.4 Non-Poissonian assumptions 

It may be more reasonable in practise to assume that the time it takes for a chunk to be 
downloaded or swapped is a random variable with a heavy-tailed distribution. This is not 
captured by our model. Indeed, the interaction times are not part of our model at all. A 
new, more detailed, model should be worked out. 

However, a crude capture of this phenomenon is the replacement of the Poisson processes 
5>£ by more general point processes, perhaps with heavy-tailed inter-event times. As long 
as these processes obey a functional law of large numbers, we can (by possibly modifying 
the scaling parameters) rephrase Theorem [1] and repeat the proof in this more general case. 
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A Drift calculation 

We consider the set of vectors 

N = {~ e F} U {e A ■ A C F} U {-e A + e A > : A C A' C F} 

U {-e A - e B + e A > + e B > : A C A' C F, B C B' C F, A' - A C B, B' - B C A}. 

For each £ £ M we define a polynomial Q^(x), by comparing ([7|) and (J2j) : 

Qe A (x) := a A 
Q- ep (x) := 5x F 

Q~e A +e A ,{x) := \ A>A '{x) 
Q-e A -e B +e A ,+e B ,(x) := fiA,B( x ) S A,A',B,B' , (13) 

where \a,A'( x ), ^a,b(x) are given by {Taj, (fT5| . respectively, and 

Sa,a',b,b> ■= HA n A' , A' - A c B, B n B' , B' - B c A). 

The variable x ranges in Z^ F ' or in M.^ F \ The algebra is the same in both cases. Define 
the drift vector field by ^2 y (y — x)q(x,y). Comparing (fT3l) and ([2]) we have 

J2(y - x )i(x,y) = ^2 (Qd x )- 

y (eM 

The latter sum appears in ([8]) , in the course of the proof of Theorem [TJ We shall verify that 

where v{x) is defined by (jl|). 

Consider the terms in the summation u(x) = X^eAf CQ(( X ) involving ( = —e A — es + 
e A' + e B'- Notice that swapping A with B or A' with B' will not change the value of 
x — e A ,x — es + e A i + ee, so we need to make sure to take into account this change only 
once in the summation. If we simultaneously swap A with B and A' with B' then neither 
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x - e A ,x - e B + e A ' + e B nor the value of Q_ eA _ eB+eAl+eB , (x) = [J.A,b(x)$a,a',b,b> will 
change because, clearly, 

VA,b(x)5a,A',B,B> = fJ-B,A( x )dB,B>,A,A>, 

as readily follows from (jlap and (|lb|) . We now see that to swap A with B without swapping 
A' with B' is impossible (unless A' = B'). Indeed, it is an easy exercise that 

5a,b,a',b' = 8 b ,a,a',b' => A' = B . 
Taking into account this, we write 

u(x) = a A eA - 5x F e F + e A + eA')^A,A'( x ) 

A,A> 

+ ^ E {-eA-e B + e A ' + e B >)fj, A ,B( x ) S A,B,A',B', (14) 
A,B,A',B' 

where the 1/2 appears because each term must be counted exactly once. The variables 
A,A',B,B' in both summations are free to move over £P n (F) (but notice that restrictions 
have effectively been pushed in the definitions of Xa,A'i I^a,b, and 6a.b,A',B')- 
Since 

E eA^A,B(x)S A ,B,A',B' = E e B fJ-A,B(x)5 A ,B,A',B', 
A,B,A',B' A,B,A',B' 

E eA'fJ-A,B(x)SA,B,A',B' = E e B > ^a,b(x)5a,b,a> ,B> , 
A,B,A',B' A,B,A',B' 

we have 

v(x) = a A e A - 5x F e F - E e A \A,A'( x ) + E e A'^A,A'{x) 

A,A' A,A' 

~ E e AVA,B(x)5A,B,A',B> + E e A > HA,b(x)5a,B,A' ,B> ■ 



A,B,A'B' A,B,A'B' 



Call the four sums appearing in this display as Ui(x),u u (x),u m (x),u lv (x), in this order. 
We use the definitions (fTa|) . (|lb|) of Xa,A'i ^a,b and find the components of the vectors 
Ui,...,u lv by hitting each one with a unit vector eg, i.e. by taking the inner products 
up = (e G , u I ),...,v^ = (e G , %}. We have: 



«?(*) = -E A ^) = -E^ G E w-G\ 1{GnA,) 

A' A' B:BdA> ' ' 

B 1 1 A' BDG 

where, in deriving the last equality we just observed that the number of sets A' that contain 
one more element than G and are contained in B is equal to \B — G\, as long as G C B: 

E 1(G □ A' c B) = |5 - G\ 1(G C 5). 
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Next, 

x A 



^E^Enr^ 1 ^) 

BDG a 1 1 



Notice that, in the last summation, G contains exactly one more element than A and is 
strictly contained in I?, so \B — A\ = \B — G\ + 1. Hence 



u ^) = P E ^Wi EA^g)^ E | g _ G i + 1 E^- g - 

1 1 A _B:_BDG 1 1 g&G 

For ti m (x), we have: 
u ?u{ x ) = ~ E ^G,b{x)Sg,A',b,b' 

B,A>,B> 

= -7E^^-E 1 ( Gc ^^- G ^)-E 1 ^^^ 5/ - ScG ) 



i ar 



-7 



B 



7 *~ xB - ( l7 ) 

As for the last term, we have: 

"ivC^) = E ^a,b{x)5 a ,g,b,b' 

A,B,B' 

= ^EE^^ 1( ^ cG ' G - AcB) 5 1(i?Ei?/ ' i?/ - 5cA) 

B A I \ I 

= 7E^TT E ^l(G^) (18) 

B geGHB 

Adding (fl"5j) and (fT7|) we obtain the first part of fl3J), while (fT7|) and (fTHj) give the second 
part. 
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